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ABSTRACT 


This study was conducted to examine the effects of discrete precipitation events on 
the short-term and seasonal evolution of ocean mixed layer temperature and salinity 
structure. This study was located at Ocean Station “P” (SO°N, 145°W) in the Northeast 
Pacific Ocean. Two numerical modeling experiments were performed. The first was to 
simulate the response of the mixed layer to hypothetical discrete (isolated) precipitation 
events. This experiment showed that the effect of a single discrete rain event can vary 
with season, with the mixed layer depth (MLD) at onset of the rain event, and with the 
strength of wind stress forcing. A single rain event can have lasting effects on mixed 
layer depth and temperature for up to 55 days after the event, depending upon the sea- 
son. 

The second experiment simulated quasi-realistic “complex” precipitation forcing, 
with a realistic distribution of synoptic events over a 13-:nonth period. For this exper- 
iment, four different precipitation “intensities” were tested. The values of temperature, 
salinity, and MLD predicted by the model were compared with observed Conductivity 
Temperature Depth measurements and with the values predicted using constant precip- 
Itation forcing. In all experiments, the modelled MLD's approximated the observed 
MLD and temperature cycle. MLD's for all cases were too deep while temperature was 
estimated well in winter and was too cool in summer. Predicted salinity provided the 
greatest discrepancy between the modeled and the observed cycle. The 96 cm, 400-day 
quasi-realistic precipitation forced case best approximated the salinity observations 
though predicted salinity was fresher than observed in winter and saltier than observed 
in summer. Model results suggest that the amount of precipitation observed at Ocean 
Station "P^ was too low to explain the observed and simulated ocean salinity and tem- 


erature Structure for the year. 
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I. INTRODUCTION 


Salinity in the oceanic mixed laver and its change in depth with time can have a 
significant impact on upper ocean dynamics. Alterations of upper ocean salinity caused 
by precipitation and evaporation through buoyancy fluxes may also change the tem- 
perature profile of the ocean, altering the sound velocity. The purpose of this study was 
to examine the effects of discrete precipitation events on the evolution of the mixed layer 
depth (MLD) at Ocean Station "P" (50*N, 145?W). 

Previous one-dimensional mixed laver prediction experiments had failed to accu- 
rately explain and predict the behavior of salinity with time and depth. Possible reasons 
for this were: 

1. the amount of precipitation and evaporation was incorrectly specified, 


2. the precipitation and evaporation distribution with time, relative to wind events 
was incorrectly specified, and or 


. the advection was incorrectly specified or neglected. 


دنا 


This study focused specifically on the variation of precipitation and evaporation with 
time. The total amounts of precipitation and evaporation are also varied. Possible ef- 


fects of advection are discussed. 


A. LITERATURE REVIEW 

Kraus and Turner (1967) provided the fundamental turbulent kinetic energy (TKE) 
theory on which many subsequent oceanic mixed laver models are based. This model 
considered the turbulent kinetic energy equation in a one-dimensional mixed laver 
scheme, using the approximately decoupled state of the mechanical and thermal energy 
equations. It was assumed that the heat input at the air-sea interface and the mass 
entrainment at the base of the mixed laver were uniform throughout the mixed laver for 
short time scales. Wind stress effects on the air-sea interface acted to increase the po- 
tential energy in the mixed layer by generating turbulent kinetic energy for mixing and 
entrainment. Entrainment was considered to be proportional to the wind energy input 
to the water column minus the work performed in overcoming the buoyant forces 
throughout the mixed laver. This type of turbulent erosion modeling considered 


buoyancy and wind forcing only, and it neglected viscous dissipation. This neglect of 


dissipation has generally been considered a serious flaw. Kraus and Turner (1967) did 
not consider salinity effects in their mixed layer model. 

Subsequent works, such as Geisler and Kraus (1969), Miropol'skiv (1970), and 
Denman (1973) have attempted to consider dissipation by assuming dissipation to be a 
fixed portion of wind stress. The Denman (1973) model was used in the first mixed layer 
simulation experiment to use Ocean Station "P^ data (Denman and Mivake, 1973). 

Miller (1976) was the first to include salinity and its effects on the density structure 
in a mixed laver model. Miller modified Denman (1973) to include salinitv and used data 
from BOMEX Period III upper ocean soundings for periods of as long as 10 days for 
forcing and initial conditions. Four specific cases were examined by selecting distinctly 
different initial conditions for the salinity and temperature profiles and forcing data from 
BOMEX surface observations. In three cases, Miller examined the effects of three dif- 
ferent initial conditions while the fourth case included a precipitation event. 

The Miller study indicated that mixed laver salinity greater than the salinity of the 
water below the mixed layer may cause convective overturning due to density instability 
unless a sufficiently large temperature decrease occurred simultaneously across the 
interface. A large density jump across the interface acted to slow the mixed layer 
deepening rate due to an increase in potential energy. Without entrainment salinity flux 
effects, the cooling induced by entrainment was greater than the heating from the sur- 
face, resulting in net mixed layer cooling. With salinity included in the model, the layer 
deepened more slowly and entrained less cold water at the bottom, allowing surface 
heating to dominate and, therefore, the mixed layer temperature to increase. Overall, the 
results of these cases demonstrated the importance of salinity fluxes and the need to in- 
clude the salinity. profile in the prediction of mixed layer deepening and resultant tem- 
perature structure. 

Miller demonstrated the potential significance of precipitation episodes in mixed 
laver modeling in a fourth case. His study demonstrated that precipitation can contrib- 
ute to salinity fluxes within the mixed layer. However, the fact that salinity regularly 
fluctuates bv +0.15% (due to horizontal advection) makes decreases in salinity due to 
rainfall difficult to detect in observations. Though Miller admittedly did not understand 
the nature of the effects of rainfall on the mixed laver, he did propose that even minor 
changes in the salinity flux due to precipitation may have a long term effect on the mixed 


laver thermal structure: 


Paulus (1978) examined the effects of a salinity profile on density structure using a 
modified version of the Camp (1976) one-dimensional mixed laver model. The effects 
of a surface salinity flux on the mixed layer were examined using a constant evaporation 
and precipitation rate. Model initialization was accomplished with hvdrocast data or 
XBT data with historical salinity profiles adjusted to realistically correspond to the 
temperature profiles. The incorporation of observed salinity profiles into the model had 
little effect on sea surface temperature (SST) and MLD during the summer regime (when 
MLD was shallow). However, when using the same forcing conditions for the winter 
regime the inclusion of a salinity profile acted to produce a warmer, shallower mixed 
laver than obtained when density was onlv a function of temperature. 

Paulus suggested that the retreat of the isothermal layer could be enhanced due to 
precipitation by slowing the deepening rate and shoaling tendency of the mixed layer 
under downward heat flux. With respect to salinity flux, Paulus concluded that imbal- 
ances in evaporation and precipitation had more pronounced effects on SST and MLD 
in the long term than in the short term. His study concluded that this was a cumulative 
effect. Therefore, Paulus concluded that the addition of salinity effects to the mixed 
laver model was desirable. 

Garwood (1977) introduced a one-dimensional bulk model of the oceanic mixed 
layer that included dynamic instability with turbulent erosion. This model modulated 
the mean entrainment rate by diurnal heating and limited maximum layer depth by en- 
hancing the effect of viscous dissipation. This allowed a cyclic steady state to occur over 
an annual time period. It recognized the important effects of long term salinity fluxes 
on the mixed layer thermal profile due to surface heat fluxes, precipitation, and 
entrainment at the base of the mixed layer. This was the model of choice for this study. 

Garzon (1987) used an oceanic mixed layer model to determine the specific effects 
of fresh water flux on the seasonal pycnocline in the mixed layer at Ocean Station “P.” 
First, sensitivity of the model proposed by Garwood (1977) was examined by performing 
a hypothetical experiment using constant forcing with a wide range of precipitation and 
evaporation values. These results were compared to an identical model run with P-E 
(precipitation minus evaporation) = 0.0. This simple hypothetical case examined the 
relative importance of excess evaporation or precipitation over a time period of one- 
year. For the case with excess precipitation, a.reduction in downward buoyancy and a 
decrease in surface salinity occurred, resulting in a shallower MLD than would be ex- 


pected if P-E=0. For the case with excess evaporation, an increase in downward 


buoyancy flux and an increase in surface salinity occurred. This yielded a deeper MLD 
than would occur with P-E=0. 

For the annual simulation at Ocean Station “P,” Garzon, using data from 1967, 
modeled temperature and salinity profiles with P-E = 0.0m'vear, P-E 7 0.5m/vear, and for 
monthlv representative values of precipitation. Overall, MLD and surface salinity from 
the model well approximated those of the actual data. However, near the end of the 
study period, MLD and salinity values demonstrated increasing error particularly in ap- 
proximating significant anomalous “spikes.” (Garzon suggested that a more realistic 
representation of discrete precipitation and evaporation events would result in more ac- 
curate model output, especially for the long term. Additionally, Garzon adjusted model 
output for advection of salinity by subtracting the linear fit from the predicted salinity 
time series. 

Precipitation and evaporation events have been demonstrated to have a significant 
impact on mixed laver depth on both short and long time scales. These factors can act 
alone or may be coupled with such events as synoptic scale storms to modify the effects 
of these events. On a short time scale, davtime heating tended to cause shallowing of 
the oceanic mixed laver due to heating of the mixed layer despite the occurrence of 
evaporation which tended to promote deepening. Upward heat flux broke down the 
mixed layer causing deepening. Precipitation during this time period promoted shal- 
lowing. thus opposing this eflect. These diurnal changes in mixed laver depth were 
moderated bv annual effects that were associated with. the annual radiation. cycle. 
Strong winds, often associated with large scale svnoptic events, when combined with low 
insolation values in winter, yielded deep mixed layer depths. However, precipitation 
events often accompanying these winds acted to lessen the effects of winds and were 
associated with discrete shallowing episodes that last for short periods of time. These 
effects were cumulative. Just as a reduction or increase in wind speed in winter vielded 
a substantially shallower or deeper mixed layer depth in summer, a reduction or increase 
of rain during the same period yielded a deeper or shallower mixed layer depth in sum- 


mer. 


B. RESEARCH PLAN: LOCATION AND PERIOD OF STUDY 
For this study, a time period was chosen to provide the best overlap of available 
data bases. The starting date of 23 November 1968 with a study period of 400-davs was 


selected for this study, since it provided data base overlap, that 1s, surface forcing and 


CTD data was available for this time period. This numerical simulation of 400-days used 
a time step of one hour. 
The study area was between 49° N and 51? N, and 144° W and 146° W centered at 


عم 


Ocean Station "P" (Fig. 1 on page 6). From 1949 to 1981 the waters around Ocean 
Station "P" have been occupied by ocean station vessels vear-round. This provided nu- 
merous data bases, such as, XBT, CTD, and atmospheric observations for analysis. 

Tabata (1961,1965) described oceanographic and meteorological conditions at Sta- 
uon "P" in detail. The dominate pressure systems located in the North Pacific Ocean 
are the low pressure svstem or Aleutian Low in winter and the high pressure system or 
North Pacific High in summer. These two systems govern the distribution of mean 
surface winds over the Northeast Pacific Ocean. In winter the Aleutian Low directs air 
flow northward into the Gulf of Alaska, while in summer the North Pacific High directs 
winds southeast at Ocean Station “P.” Additionally, monthly mean wind speeds are 
approximately twice as large in magnitude (13m/s) 1n winter as in summer (6m s). 

Station "P" is located in a distinct oceanic regime (Fig. 2 on page 7). The surface 
waters of the region are under the influence of the weak westward flow. Specifically. this 
flow is composed of the West Wind Drift and the recirculated waters from the Alaska 
Stream. The former carries mixed waters of the Kuroshio and the Ovashio, which are 
generally indistinguishable by the time that they reach Ocean Station “P.” To the east 
of this area, the current turns north, forming the Alaskan current. 

The vertical structure of this water 1s characterized by three zones: (1) the upper 
“seasonal” zone (0-100m), (2) the halocline (100-200m), and (3) the lower zone (> 200m). 
Within the upper zone the salinity is relatively low (~ 32.7ppm). With respect to sea- 
sonal effects, in winter the temperature and salinity in zone | are relatively homogeneous 
and isothermal. In summer a significant vertical salinity gradient exists in this upper- 
most zone, as Well as a well defined thermocline, with a temperature decrease with depth 
of as much as 8°C per 20m. In the halocline, zone 2, salinity increases with depth by 
as much as 1.0ppm per 100m. In zone 2 temperature decreases with depth typically, 
unless an inversion is present. In the lower zone salinity increases only gradually with 
depth, while temperature decreases gradually with depth (Fig. 3 on page 8). 

Tabata (1961) indicates that there is an annual average excess of precipitation over 
evaporation of ~ 0.5m. Garzon (1987) also estimates an excess annual precipitation. 
Despite this annual excess of precipitation, salinity within the water column does not 


decrease on an annual time scale, although it can vary as noted by Tabata (1961), 


Region Surrounding Ocean Station "P" (from Tabata, 1961) 
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Fig. 3. Example Temperature and Salinity Profiles 
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especially within the more dynamic mixed layer. Horizontal advection must account for 
this readjustment. As noted by Tabata, upwelling is weak and nearly insignificant within 
these waters, when compared with horizontal advection. As a result this salinity balance 
appears to be maintained by horizontal advection as more saline water is transported 
into the area by ocean currents. 

Despite the consistent westerly current coming into this area, anomalous events 
have occurred resulting in horizontal transport of water into the area from a non- 
westerly direction. Tabata (1961) noted an unusual increase in temperature and decrease 
in salinity. within the upper zone and halocline during spring 1957 through summer 1958. 
This was attributed to transport in water lving south of Ocean Station "P" into the studv 
area. The suggested cause for this transport was the shift of winds in the region associ- 
ated with unusual location of the atmospheric pressure svstem over the northeast Pacific 
Ocean. This displacement of the North Pacific High vielded westerly instead of the usual 
southeasterly winds, which could have generated southerly drift currents in the vicinity 
of the station. 

Overall, the currenis in this area are weak. Except for occasional anomalous trans- 
port associated with unusual synoptic weather svstems, salinity follows a roughly annual 
cvcle maintained bv horizontal advection, and the waters of the study area are relatively 
homogeneous horizontally. These factors combine to provide an excellent region for a 


one-dimensional modeling study of the thermodynamics of the upper ocean. 


ll. DATA SOURCES AND PREPARATION 


A. DATA SOURCES 

This study was performed using various sources. Atmospheric forcing data was or- 
iginallv obtained from the National Weather Records Center, Asheville, North Carolina. 
This data included threc-hourly observations of wind direction and speed, air temper- 
ature, sea surface temperature, dew point temperature and cloud coverage in octaves. 
Atmospheric data was interpolated from three hourly to one hourly data for computa- 
tion of hourly surface fluxes. Precipitation forcing data was not directly available and 
was derived as explained later in this chapter. 

Two types of oceanic data were used. CTD/STD data was provided by the Marine 
Environmental Data Service, Ottawa, Canada. These salinity and temperature casts 
were digitized from depths of 0 to 1500 m at various intervals, with more concentrated 
data provided within the mixed layer. These casts were not taken at regular intervals, 
but were taken as operationally feasible. There were 119 casts available for analysis 
during the study period from station “P.” Bathythermograph casts were not available 
for this time period. 

Three hourly atmospheric forcing data was used in this study to compute the fluxes 
of heat, radiation and momentum necessary for forcing the mixed layer model with an 
hourly time step. Specifically, these boundary conditions were wind stress in two- 
dimensional components, incident solar radiation, heat loss (long and short wave), and 
precipitation and evaporation. The calculation of precipitation will be discussed later. 
Temperature and salinity profiles from STD taken 23 November 1969 were used to ini- 
tialize the model profiles, and CTD. STD casts were used for verification. Bucket tem- 


peratures were used for comparison of with predicted temperature. 


B. VERTICAL INTERPOLATION 

The oceanic profiles, from both CTD/STD sources were vertically interpolated from 
O0 to 199 m to a standard grid of 1 m spacing. Mixed laver depth, mixed laver temper- 
ature, and mixed layer salinity were derived from these interpolated profiles. There is 
not a standard way to compute mixed laver depth from a CTD/STD profile. For this 
study the mixed laver depth has been defined to be that depth at which the temperature 


profile changes exceed 0.10* C in 5 m. Mixed laver temperature and mixed laver salinity 
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values were obtained by averaging temperature from 1 m to the mixed laver depth as 


determined above. 
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lil. THEORY 


A. OCEANIC MIXED LAYER STRUCTURE 

A familiarization with the definitions associated with the oceanic mixed layer ther- 
mal structure 1s required to understand the dvnamics involved in mixed layer modeling. 
As defined in Garzon (1987) the oceanic mixed boundary laver 1s assumed to be the fully 
turbulent region of the upper ocean bounded above by the air-sea interface and below 
bv a relatively static water mass. The mixed layer is considered to be “quasi- 
homogeneous," that is, it is approximated to be well mixed. Thus temperature within 
the mixed laver is considered to be relatively constant. Fig. 4 on page 13 graphicallv 
depicts an oceanic vertical temperature profile. The turbulent mixed laver of depth h is 
contained within the upper most oceanic region, below which lies the turbulent 
entrainment zone, from depths z=-h to z=-h-d, with z positive upward. Within the 
entrainment zone there is a temperature jump of AJ. Below the entrainment zone is the 
quasi-stable water mass and here temperature decreases with increasing depth. 

Mixed laver structure (Fig. 5 on page 14) for salinitv is similar to that for temper- 
ature (Fig. 4 on page 13). There is a relatively constant salinity within the turbulent 
nuxed laver from z=0 to z=-h, a salinity jump AS across the entrainment zone from 

=-h to z=-h-6, and a gradually increasing salinity with depth below the mixed layer. 

An oceanic mixed layer prediction scheme requires an accurate (or near accurate) 
assessment of the deepening or shallowing rate of that laver. Garwood (1977) accounted 
for both of these processes. The equations used by Garwood to best describe these 


processes are developed as follows. 


B. ENTRAINMENT 
The Garwood (1977) one-dimensional mixed layer model considers both turbulent 
erosion and dynamic instability. Fundamental to the solution of this model is an accu- 


rate assessment of the entrainment velocity, 1۲.۰ Starting with: 
Ch j ; 3 
a ٥٢ (3.1) 


Assuming no mean vertical motion, W (-h) 7 0, this equation becomes: 





oh 
W,-- ! 
a (3.2) 





Fig. 4. Example Temperature Profile 


The solution for IV, can be found by solving the total turbulent kinetic energy (TKE) 


budget written as: 
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where: 
E=U'"+V"+ Wa is the Total TKE. 
(see Table 1 on page 20 for definitions of constants) 
The five terms in Equation (3.3) from left to right are: 


l. the time rate of change of turbulent kinetic energy, 





Fig. 5. Example Salinity Profile 


2. shear production, 
3. buovancy flux, 
4. turbulent diffusion (or redistribution), and 


5. viscous dissipation. 


This equation may be examined for the entrainment zone. Within the entrainment zone, 
dissipation and shear production are assumed negligible. This produces a mean turbu- 


lent kinetic energy equation for the entrainment zone of the form 
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across the entrainment zone of thickness ó. Integrating equation (3.4) vields: 


د 





ele 
to 


=h zh جج‎ 
| 2 تو‎ | xgT'Waz — | n E. | (3.5) 
Z Po lS 
h-ó —h-ó 


Assuming steady state, the time rate of change of turbulent energy may be neglected. 


Equation (3.5) becomes: 


—h / —h 
ve - vt E P -- ^ 
| ern dz | H | oa ji. 7o (3.6) 


Solving this integral produces equation (3.7). 


DN تس شر‎ E : 
— (TH | ¿+ FIL o) 1" ۷ال مہ + ہے‎ (3.7) 


When specifically examining this equation at the base of the mixed laver (or the top of 
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the entrainment zone): 


P" 


sgT TW = E mó 


p 
nb 3.8 
po ) | —h ( ) 
Kraus and Turner simplified the solution at this point by neglecting the redistribution 
term and combining the entrainment mean kinetic energy equation with the steady state 
turbulent Kinetic energv equation for the mixed layer. Garwood does not neglect this 


term, but estimates its magnitude with the following parameterization: 


6. See. OP > جام‎ 
A A = = p 

Gz ( 2 Po Ma = (3.9) 
Thus the redistribution term is equal to the mean turbulent energy available near the 
base of the mixed layer divided by the time scale required to transport this energy into 


the area. Parameterizing this time scale as 


h 


eee os 3.10 
| H7 | ( ) 


Te 
where | W'| is the root-iean square vertical turbulent velocity, | W/| 2 (W2)'2, and the 


time scale is approximated as the absolute value of the mixed layer depth divided bv the 


absolute value of the vertical velocitv. 
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Using Equations (3.8 to 3.10) and assuming that T'W'| , 2 ATW, , 


| 11 ]| < ل ب 


gluAT — BAS) W, = 7 


En. 
Here vertical transport of turbulence exactly balances entrainment buovancy flux. This 


equation may be solved for the entrainment velocity: 


/! F 
"۲ IW|<E> 


* — hg(a AT — BAS) Ca 


To solve this equation for W, the values of < E >, | W'| , T, h, and AT must be 
known while g and o are constants (Table 1 on page 20). From Garwood (1977) the 


total turbulent Kinetic energy is solved by using: 


D |o» 


(h< E >)= G + 8ر8 -و8 - ,م0‎ (3.13) 


In order to solve equation (3.13), the following equations for shear production at 
Z=0 and z=-h, and buovancy damping’production at z=0 and z=-h, and dissipation 
are: 


wind shear production, 


m t 7 
C= ادا‎ y ms (+ p" (3.14) 
entrainment shear production, 


G= [< E> +< i >] (3.15) 


where the x and v components of the mixed layer current are: 
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buovant damping, production due to surface buovancy, 


Qo c 
zd — pgS(E — P), ono) 


buovant damping due to entrainment, 
B, = oghAT W,, and (s) 
dissipation, 
D=2< E>”. (3.18) 


C. MIXED LAYER THERMODYNAMICS: TEMPERATURE AND SALINITY 
BUDGETS 

Wind stress generated turbulence at the air-sea interface provides the energv to de- 
stabilize and entrain the water below. This mixing process deepens the mixed laver and 
mixes temperature and salinity between the surface and the entrainment zone. Thus, the 


conservation of heat requires that: 





cr FP ae” (3.19) 
Gt 02 
Vertically integrating Equation (3.19) across the entrainment zone, z=—hÁh-—ó to 
z — — h, and assuming 7" FV” vanishes below z = — Ah — ô gives 
— AUAM. — UE — hy) (3.20) 
where: 


r= T(z =—h)— T(z =—h—8). 


Vertically integrating Equation (3.19) across the entire mixed layer, z=-h to z=0 


gives 
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Using Equation (3.20) and assuming T to be well mixed in z, Equation (3.21) becomes 
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Here: 
۰٠۹77 و‎ +١ ٤ 
2 ات‎ =- T T W'(0). 
The one-dimensional salinity budget for the oceanic mixed laver is: 
د‎ (3.24) 
Ol OZ 


It 1s assumed that: 


]. S(z) 1s well mixed from the bottom of the mixed laver to the top of the mixed laver, 


ty 


S(z) has a near discontinuity between z= -h and z= -h-ô, and 


3. the vertical flux of S (or S' I^ ) vanishes at z < —A — Û. 


Fig. 5 on page l4 depicts this condition. Here the turbulent region within the mixed 
layer yields a near homogeneous state with respect to salinity. Examining Equation 


3.24) and integrating gives: 


© 





Ud 


AOS چا‎ 
> dz = اک‎ 2 dz (3.25) 
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With h variable in time, application of the Leibnitz Rule of integration and assuming ۵ 


is vanishing small, 


AS 


ات 


=P S- HFE S (3.26) 


Assuming that turbulence is negligible below the mixed laver, ie. W'S'( —h— 6) —0 


This makes the “jump condition” for salinity entrainment into the mixed layer: 


As == — W"S'( SR). (3.27) 


where: 
AS 2 S(z2—h)—S(z-—Rh-—6) 


D. SALINITY EFFECT ON MIXED LAYER DYNAMICS 
l. The Shallowing Effect 
Kraus and Turner's (1967) model provides an approximation for the non- 


entraining or shallowing mixed layer, W, =0, 


6 U} 


D 3.28‏ اچچ تو را اک 
C 0800/(0 e‏ 


Here L is the Obokov length scale and equation (3.28) 1s the solution when there 1s no 


entrainment. The surface buoyancy flux is B, = (agq,/(pC,)) and U. is the surface friction 


velocity, with U? —(p,Jp)C;U? . The value of h is inversely proportional to Q,. The 
mixed layer temperature is even more sensitive to Q, than is h: 
CT 2 1 4 
"ar 7 Qu(o Ch) oc Q' oc p (3.29) 
1 


If the buoyancy flux attributable to precipitation and evaporation 1s included, Equation 
(3.28) becomes 
we E (3.30) 
fü DE Gl mar E و‎ $9 
Cy agQgl(pC,) + Bg(P — E)S 
Where the salinity concentration coefficient, f = 0.76x10-*Kg/gm , and P-E is the net 
precipitation minus evaporation in m,sec. With the salinity flux included, the surface 


buoyancy flux becomes: 
By = ag Qy PC, + Bg(P — E)S, (3.31) 


and the effect on mixed layer temperature becomes: 


6 

ar و18020)02/ |( م/و0) ح‎ + Bia(P — EJS). (3.32) 
As a result, precipitation minus evaporation influences mixed layer temperature indi- 
rectly. For positive Q, (net heat into the surface), a positive (P-E) will decrease h, con- 
centrating the heat into a shallow layer. For the case where evaporation exceeds 


precipitation a decrease of mixed layer temperature with respect to time results. 
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2. The Deepening Effect 
For entraining, deepening mixed lavers, the buoyancy discontinuity at the bot- 


tom of the mixed laver must be augmented to include a possible salinitv jump: 
Ab = agAT — ۵5ع‎ (3.38) 
WV ene: 
AS =S —S(-h-5), 
S is the mixed layer salinity, and 
S( — h — ó) 1s the salinity just below the mixed layer. 
The mixed laver salinitv budget 1s: 


es | S(E- P) | W,AS 
= چک‎ (3.34) 
CI h h 
Equation (3.34) demonstrates that the change in salinity with respect to time 
depends on fresh water flux at the surface, entrainment velocity, salinity advection, 


mixed laver depth, and upon the salinity jump at the base of the mixed laver. 


Table 1. MODEL CONSTANTS 


Model Constants 


a thermal expansion coefficient 
salinity ات ۷۳ چ‎ E ~Š x 10 * ppm! 


[tuning coefficient | = E 
[tuning coefficient — | — 3 
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IV. HYPOTHETICAL EXPERIMENTS 


A. EFFECTS OF DISCRETE PRECIPITATION EVENTS 

The purpose of this study was to determine the lasting effects of discrete precipi- 
tation events on the mixed layer. For “case 1” in this study short-term precipitation 
events of 10 cm in one day with no evaporation was applied monthly to the Garwood 
(1977) model over a 365 day study period. This vielded twelve separate events for de- 
tailed study. At all other times during the year both evaporation and precipitation were 
assumed to zero. The resulting output from the model was examined individually for 
30 day periods following each rain event to better reveal long term trends. The precipi- 
tation event was applied on the fourth day of each monthly period. Actual temperature 
and salinitv profiles taken from a CTD cast for November 23, 1968 were used as initial 
conditions. Heat and momentum fluxes were computed from the actual surface 
meteorological and observations. 

Fig. 6 on page 22 and Fig. 7 on page 23 are a conglomeration of the results of these 
twelve separate events. Other figures in this chapter represent the individual events that 
were selected for closer examination at a 60 day period. 

The effects of precipitation events on the MLD are graphically represented by the 
figures in this chapter. The equations developed in Chapter 3 will be used to explain the 


results and are repeated here: 
09060: 5 (4.1) 


The change in mixed layer temperature (T) with respect to time is a function of: 
EM LID h, 
2. entrainment rate, iV, 
3. temperature jump at the base of the mixed layer, AT, and 
4 


. net downward heat flux, و0‎ 


For a case of no entrainment, WM”, =0, the temperature will rise for Q, > Q, and it will fall 
for Q, «0. The rate of rise and fall is inversely proportional to h. If Q, x 0 and there 


is entrainment, IF, > 0, the mixed laver temperature can only decrease if AT > 0. The 
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Fig. 7.  Conglomeration of Six Events (Jun to Nov 1969): (a) diflerential MLD, 
(b) diflerential temperature, (c) differential salinity, (d) MLD with E-P=0, 


(e) wind speed, and (f) net surface heat flux. 
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rate of decrease is greatest for small h. Various factors effect IW., and will be considered 
later. 

Two model predictions may be compared. In this experiment the first case included 
the discrete precipitation event, while the second case had no events, with E-P =0 for the 
entire model integration. When subtracting the results of one study using equation (4.1) 
from another such study for the same time using the same equation, the results provide: 

1. differential MLD (A, — h), 
2. differential temperature (T, — 7,), and 


3. differential salinity (S, — S,). 


Differentiating the differential temperature with respect to time yields: 


ATT) (1 130 [ATI ATM 4 
Ol 1 hy, h, ) Des hy hy 


Therefore, in order to provide a value other than zero for ¢(T, — T,)/ér , the differential 


temperature, one or more of the following conditions must exist: 
Lose 
2. MÁ, 5 Wes, and 
DUNT AT. 


Again, considering the many variables that may alter the entrainment rate, there are 
numerous factors which may effect this differential temperature, making interpretation 
of even a simple precipitation event complex, but more easily examined than the case 
with a more realistic and complex E-P time series of forcing. 

Initially, in all twelve events with all forcing values equal (for the first five davs), the 
differential temperature ôT = T, — T, must equal zero. This can be seen in Fig. 6 on 
page 22 and Fig. 7 on page 23 for the first five days of each 30-day period. Upon initi- 
ation of the precipitation event at day five of each events A, = A; and AT, 2 AT; , but 
W, 1s instantaneously altered in the precipitation case. In this case with E-P < 0,there 


results in an immediate increase of surface buoyancy, B, , 





Q "s 
By = ag L — figS(E — P). 
7" 


An increase in D, acts to decrease the total TKE from equation (3.13), 


e 
p, (h< E>)= G+ و8 - ريا‎ - B, — D 


and, therefore, entrainment rate from equation (3.12), 
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As W, decreases, so will MLD, which acts dynamically to shallow the MLD relative 
to the case with P-E=0. Thus a precipitation event increases the downward surface 
buoyancy through decreasing salinity, which ultimately results in a shallower MLD. 
This initial result should be apparent throughout the twelve experiments where there 1s 
precipitation. 

While the immediate change in MLD through reduced entrainment 1s expected on 
a short time scale (during the rain event), longer term changes in mixed layer dynamics 
may result from a single short-term precipitation event. In response to mixed layer 
shallowing, MLD and the temperature jump at the base of the mixed layer will change 
from that of the non-precipitation case. At this point all three possible conditions re- 
quired to yield an increase in the differential temperature may exist. 

Thus analvsis of the effects of a discrete precipitation event requires a case by case 
examination. Anv differences in MLD will vield a non-zero magnitude of |/ — 7l 
which will effect the differential temperature. A fluctuation in differential MLD will also 
moderate the effects of Q, , AT, — AT;, and W,, — ÀW,,. Any difference in the temper- 
ature jumps, |AT, —AT,| #0 , will also yield an associated change in the differential 


temperature and a transient response, as well as the altered value of IV, . 


B. ANALYSIS OF PRECIPITATION EVENT EFFECTS 
I. Event #1 (Dec 1968) 

Fig. 6 on page 22 represents the MLD output for the model with E-P=0. This 
figure indicates that this event began in the late-fall season when the mixed layer was 
at 80 m and deepening. Fig. 6 on page 22 and Fig. 8 on page 26 indicates that there 
was a net heat loss (Q, < 0) during this period associated with a mid-latitude, winter re- 
gime. There were moderate winds during this time period (U = 5.0m/s) , and maximum 
values of U,,, = 20.0m/s which occurred immediately before, during and after initiation 


of the precipitation event of day 5. 
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Fig. 8. Event #1 (Dec 1968): (a) differential MLD magnified from day 2 to day 


6, (b) difTerential temperature, (c) difTerential salinity, and (d) wind speed 


Prior to initiation of the precipitation event all differential values equaled zero 
(Fig. 6 on page 22 and Fig. 7 on page 23). This was true for all twelve events. Upon 
onset of precipitation at dav four, MLD decreased from 65m to 10m, vielding a decrease 
in the differential MLD. This was coincident with shallowing of the mixed laver asso- 
clated with a lull in wind. Approximately 12 hours after precipitation onset, however a 
moderate wind event followed immediately by another of magnitude U~ 5.0m/s oc- 
curred. These increases in wind stress occurred during the precipitation event and re- 
sulted in two increases of the differential MLD and temperature (points a and b). Thus, 
the layer deepened, due to decreased W, causing the layer to cool. 

Immediately following events a and b the wind decreased to near zero causing 
a large decrease in MLD. Here the differential MLD reached a maximum, as the MLD 
in case 1 shallowed to near the surface. Associated with this decrease in MLD in case 
l, was a negative differential temperature. This was the result of the shallower MLD in 
case | allowing for the increased heat loss bv the layer with Q, >Q (net heat loss) in 
winter (point c). The differential MLD increased then decreased, as induced by the 
onset and drop-off of a wind event (points d and e). Again this produced an associated 
positive and negative differential temperature, respectively. At day six a strong storm 
occurred, (U ~ 13.0/s) (event f) followed immediately by a series of larger fluctuations 
in wind speed. These storms brought the MLD of case 1 to near that of case 2. In- 
creased buoyancy, causing decreased entrainment, caused the MLD in case | to deepen 
less than that of case 2 accounting for this significant differential temperature event at 
day six as case 1. This entrained less cold water than in case 2. After day six a series 
of winds (U ~ 5.0m/s) continued to mix the layer, bringing the differential MLD’s and 
temperatures closer to zero. Although the effect of the precipitation event on the mixed 
layer decreased with time, the layer temperature and depth were effected measureably for 
this experiment period of 60 days. 

There was an immediate decrease in salinity as precipitation diluted the layer, 
and this decrease continued until cessation of the precipitation event. Mixing by 
entrainment decreased the differential salinity to a value of -0.4ppm, after which time it 
rapidly approached that of case 2, aided by wind mixing. 

2. Event #2 (Jan 1969) 
As in Dec 1968, prior to the precipitation event the MLD was relatively deep 


(h ~ 90m) and continued to deepen during this winter regime. The values for Q, still 


27 


indicated a net loss to the atmosphere and strong winds were predominate 
(U ~ 5.0m/m) . 

After the onset of the precipitation event, the results of the event were very 
similar to those of event 41. Therefore, no specific figures were provided other than 
those found between days 30 and 60 in Fig. 6 on page 22. 

At the beginning of the precipitation, differential MLD generated a large de- 
crease ( ~ — 100.0m) as the MLD of case 1 nearly reached the ocean surface. This was 
concurrent with a wind speed lull for two davs, enhancing this MLD shallowing. Asso- 
ciated with this MLD decrease was a decrease in the differential temperature as the 
shallower MLD in case 1 enhanced the temperature decrease associated with a net loss 
of heat in winter, Q, « 0. By day 40 a series of strong wind events (U ~ 7.5m/s) acted 
to converge the MLD of case 1 back to near that of case 2. Likewise, the differential 
temperature approached zero in a stepwise fashion as the temperature in case l ap- 
proached that of case 2. Decreased entrainment rate due to increased buovancv did not 
allow the MLD of case 1 to remain as deep as that of case 2. This caused the temper- 
ature of case 1 at dav 25 to exceed that of case 2. It did not, however, exceed the in- 
crease seen in Event #1 as the difference in the MLD's between case ] and case 2 were 
smaller in this event. : 

There was an immediate decrease in salinity as precipitation diluted the layer 
and this decrease continued until the end of the rain event. Mixing by diffusion and 
advection decreased the differential salinitv to a value of -0.3ppm, after which time it 
rapidly approached that of case 2. 

3. Event #3 (Feb 1969) 

The MLD prior to the precipitation event was deep (/ ~ 95m) and deepening 
during the winter regime. The value of Q, still indicated a net loss of heat to the at- 
mosphere, but the magnitudes of Q, and of the wind speeds were smaller than for the 
previous two events. 

After onset of the precipitation event in case l, the results of the event were very 
similar to those of Event #1 and #2. As in Event #2, no specific figures were provided 
other than those found between days 60 and 90 in Fig. 6 on page 22. As the rain began, 
differential MLD generated a large negative event (~ — 70.0) as the MLD in case 1 
nearly shallowed to the ocean surface. This was concurrent with a wind lull of four days, 
enhancing this shallowing. Associated with this MLD decrease was a decrease in the 


differential temperature as the shallower MLD in case 1 increased the temperature drop 
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associated with a net surface heat loss. During the end of the precipitation event there 
was a small peak in both the differential MLD's and temperatures associated with a wind 
speed maximum. By day 70 series of strong wind events (U ~ 5.0m/s) caused the MLD 
of case 1 to converge to near that of case 2. Likewise, the differential temperature ap- 
proached zero in a stepwise fashion as the temperature in case 1 approached that of case 
2. Decreased entrainment rate due to increased buovancy did not allow the MLD of 
case 1 to remain as deep as that of case 2. This caused temperature of case 1 at day 70 
to exceed that of case 2. It did not, however, exceed the increase seen in Events #1 and 
Z2, as the difference in the MLD's between case 1 and case 2 were smaller in this event. 

There was an immediate decrease in salinity in the first case, as in the first two 
events, as the precipitation diluted the layer. This decrease continued until cessation 
of the precipitation event. Mixing by entrainment decreased the differential salinity to 
a value of -0.3ppm, after which it rapidly approached that of case 2. 

4. Event #4 (Mar 1969) 

Examining the MLD model output for the E-P=0 case (Fig. 6 on page 22), it 
appears that this was the period of deepest MLD (A > 100m). This was also the begin- 
ning of the period of transition from the colder winter regime to the warmer summer 
regime. Here the MLD began to shallow due to decreasing wind speed and Q, > 0. 

This event provided different results than Events #1 to #3, and Fig. 9 on page 
30 is presented for better resolution. As in the previous events, there was an initial de- 
crease in both the differential MLD and in temperature (point a). This was associated 
with the precipitation event and a large decrease in U from 18.0 to near 0.0 m's. During 
this shallowing period there was a slight wind increase causing U to peak to 5.0 m/s, and 
a deepening in both case 1 and case 2. This was less in case 1 than in case 2. This dif- 
ferential in MLD generated a positive differential temperature (point b). The positive 
differential temperature immediately decreased as the wind dropped to zero, and the 
MLD became very small, order 1 m (point c). Net heat loss cooled the upper layer in 
case ] more than in case 2, generating a severe negative spike in the differential tem- 
perature (^ — 0.10*C). A series of mild wind speeds until day 98 yielded various negative 
and positive peaks in the differential MLD's and temperatures. 

At day 100 a strong wind (U ~ 8.0m/s) returned the MLD in case | equal to that 
of case 2, making the differential MLD zero (point d). This mixing decreased the dif- 
ferential temperature to near zero, but it remained slightly colder as the entrained water 


from below the mixed layer remained warmer than water of the previous shallow laver. 
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Fig. 9. Event ff4 (Mar 1969): (a) diflerential MLD, (b) dillerential temperature, 
(c) differential salinity, and (d) wind speed 
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Strong winds at days 113, 130, and 140 continued entrainment, thus, diffusing the colder 
layer in case | with the warmer water below, and bringing the differential temperature 
closer to zero (points e to g). 

Near the end o^ this 60-dav experiment there was a very strong mixing event at 
day 147, induced by a U peak of 15.0 m/s with a duration of 3 davs. Almost 60 days 
after the single precipitation event, there was still an effect on the mixed layer. Exam- 
ining the differential MLD and temperature, 1t appeared that the MLD in case 1 did not 
deepen as much as in case 2. This deeper MLD caused the temperature in case | to begin 
to exceed that of case 2 (point h). 

Following the cessation of precipitation, in the first case, diluted salinity was 
mixed by entrainment, decreasing the differential salinity to a value of -0.4ppm. After 
this it rapidly approached that of case 2 by wind mixing. 

5. Event #5 (Apr 1969) 

Near the end of the March studv period (Event #4) the winds decreased and the 
MLD shallowed (Fig. 6 on page 22). After the wind speed lull from day 123 to dav 127, 
there was a steady rise in wind to the peak value at day 147 of U~ 14.0m/s . This 
gradual increase provided the TKE to mix the laver to to the deepest point for the year 
(h^ 125.0m). This was followed bv a period of shallowing, induced by an abrupt fall in 
wind speed (Fig. 10 on page 32). The value of Q, was positive and increasing. net heat 
flux. This served to provided a net heat gain in this summer regime, as opposed to the 
heat loss that occurred with the winter regime. 

Fig. 10 on page 32 demonstrated the initial decrease of MLD in case l, to be 
greater than that of case 2 (point a). Due to the low wind speeds coincident with the 
initiation of the precipitation event, the mixed layer of case 2 was also relatively shallow 
(A ~ 27.01) . Thus the increase in the differential temperature was less than that of one 
day later (at point b) following a wind-driven deepening event. This produced a large 
decrease in the differential MLD. A series of shallowing and deepening events followed 
and were associated with differential temperature adjustments. This brought the differ- 
ential MLD and the differential temperature back toward Zero. 

The largest decrease in differential salinity in Fig. 6 on page 22 occurred in this 
event. In this event the precipitation event was initiated after a long period of low 
winds, providing a shallow MLD prior to onset of the precipitation event. This may 
have been an unrealistic example. Given that increased wind speeds generally ac- 


company large precipitation events like that of day 124, this hypothetical event would 
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Fig. 10. Event #5 (Apr 1969): (a) diflerential MLD, (b) differential temperature, 


(c) differential salinity, and (d) wind specd 


have been unlikely to occur, but it permits an interesung examination of model dvnam- 
ICS. 
6. Event #6 (May 1969) 

Between day 150 and dav 180 the MLD in case 2 was similar to that of Event 
Z5, but more shallow. This shallowing was the result of decreased wind forcing as seen 
in Fig. 6 on page 22. Met heat gain to the ocean surface continued to increase as sum- 
mer approached, and some of the largest net downward heat fluxes occurred during this 
per10d (Ob max 70017). 

The dynamics of this event resembled those of Event #5. No 60-day period fig- 
ures Were provided. Fig. 6 on page 22 provides adequate resolution. As before, there 
was an initial shallowing producing a positive differential temperature associated with 
the shallower MLD and the summer heating regime. Until day 165, the differential 
MLD and the differential temperature resembled those of the previous events. After day 
165, however, the resulting differential temperature profile differed from those of the 
previous events. Here, the differential temperature reached a post-precipitation mini- 
mum value of 0.05°C and continued to increase slowly through the rest of the study 
period. 

This steady increase in differential temperature, when compared with the tem- 
perature plot of the mixed layer, indictated that the mixed layer was warming faster for 
case 1 than case 2. The MLD's were equal in depth during this period, so that this 
condition could not have been due to case | being shallower than case 2. The possible 
causes Were decreased H',, due to decreased salinity, and/or increased AT. An increased 
AI could be caused by the mixed layer in case 1 being warmer than that of case 2, and 
the temperature of the water mass below being the same for both cases. If so, in the 
summer regime this steady increase in the differential temperature would appear to have 
been “self-perpetuating.” 

There was an immediate decrease in salinity in the first case, as in the other 
events, as the precipitation diluted the layer. This decrease continued until cessation 
of the precipitation event at day 156. Mixing by entrainment decreased the differential 
salinity to a value of -0.3ppm, after which time it rapidly approached that of case 2. 

7. Event #7 (Jun 1969) 

Though the MLD for case 2 began relatively deep (h ~ 70m) it decreased to the 

shallowest mean period value (h ~ 20m). This was caused by an initial wind event at day 


181 (U ~ 6m/s) , and followed by a relative lull from day 182 to 205. Another wind 
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Fig. 11. Event 47 (Jun 1968): (a) differential MLD, (b) differential temperature, 
(c) differential salinity, and (d) wind speed 


increase occurred, similar in magnitude to the first at day 205 (Fig. 11 on page 34). The 
net heat flux remained stronglv downward, but not as strong as in the previous study 
llus. 7 on page 23). 

Initially, the MLD for case 2 was shallow (/5 — 11m) . At the beginning of 
precipitation, the MLD of case 1 shallowed to only 3m, yet there remained a small dif- 
ferential MLD, producing a small differential temperature « 0.1^C (Fig. 11 on page 
34,point a). This difference decreased with continuing low level winds which acted to 
deepen the MLD in both cases and vielded a decrease in the differential temperature to 
near zero (point b). 

The most significant deepening event of this period then occurred (point c), 
which produced the largest differential temperature of the period ( - 0.5%C). The mixed 
laver of the rain case then behaved as expected, stepping differential MLD and differ- 
ential temperature back toward zero. The differential temperature for this period was 
much larger than for the other periods studied. This appeared to be due to the very 
Shallow MLD’s achieved by case 1, enhancing the effect of the downward summer heat 
flux. 

The shallow MLD at the onset of precipitation did vield a stronger than normal 
differential salinity (~ —0.070ppm). The differential salinity acted as in the previous 
events, returning to an equilibrium value of ~ — 0.02ppm . 

8. Event #8 (Jul 1969) 

This event coincided with a shallow (h ~ 30m) trend in the MLD for case 2. 
This was associated with low wind speeds (U < 5m/s) (Fig. 7 on page 23 and Fig. 12 on 
page 36). Net heat flux appeared to be steady (Q, ~ 1801) . 

There was an initial positive increase in differential temperature which occurred 
in a stepwise manner, commensurate with the differential MLD (Fig. 12 on page 36). 
This occurred. until wind speed decreased to zero at day 216 (point a). From day 216 
to 226 (points a to b) the differential temperature remained nearly constant at 
^. 0.065*C, with only a slight increase in time. The MLD in case 1 equaled that of case 
2 (and the differential MLD ~ 0) in this period. This condition was simular to the slight 
increase which was note in for event #6, and the physical explanation is the same. 

After day 226 (point b) there was a large increase in the differential MLD and 
differential temperature cue to a U spike of ~ 4.3m/s. This acted to cool the mixed layer 
of case 2 more than that of case 1 and yielded a positive increase in the differential 


temperature (~ 0.1°C). The differential temperature then fluctuated around 0.22?C for 
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Fig. 12. Event #8 (Jul 1968): (a) differential MLD, (b) differential temperature, 


(c) dilTerential salinity, and (d) wind speed 
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the remainder of this 60-day period as mild winds allowed shallowing of the mixed layers 
in cases ] and 2 with intermittent weak spikes in wind speed producing some oscillation 
around this mean. 

There was an immediate decrease in salinity in the first case, as in the other 
events, as the precipitation diluted the laver and this decrease continued until cessation 
of the precipitation event. Mixing by entrainment decreased the differential salinity to 
a value of -0.3ppm, after which time it rapidly approached that of case 2. 

9. Event #9 (Aug 1969) 

With the exception of the initial small increase in differential temperature, this 
event produced results very much like that of Event #8. This small fluctuation was ac- 
counted for bv the weak wind mixing (U ~ 4y/s) occurring upon onset of precipitation. 

10. Event #10 (Sep 1969) 

During this study period the MLD in case 2 was relatively shallow, and begin- 
ning to deepen, as there was a seasonal shift from summer to winter (Fig. 7 on page 
23). This was associated with steadily increasing wind speeds (Fig. 13 on page 38). As 
the seasonal regime began to shift to a winter regime, the Q, shifted from one of net 
downward heating toward one of cooling, although, during this period net heating re- 
mained the slightly positive. 

As with the results of Event #9, these experimental results resembled those of 
Event #8 without the earlv fluctuating. In this event, prior to precipitation onset, the 
MLD in both cases had shallowed to near surface due to a wind speed lull from dav 273 
to 275 (Fig. 13 on page 38). A slight fluctuation at dav 274 (point a) generated a small 
downward event in differential MLD (~ — Im) which promoted a slight upward differ- 
ential temperature. With a follow on U event at dav 275 (point b) of ~ 4m/s, the largest 
of all positive differential temperature in the twelve events was generated. This occurred 
along with a differential MLD of ~—25m. Wind events followed which increasingly 
deepened the MLD of case | and brought the differential temperature back to near zero, 
and constantly decreasing. 

This event produced the largest differential salinity of all, concurrent with this 
anomalous differential temperature fluctuation. Again, the very shallow MLD estab- 
lished before precipitation onset appeared to have contributed to this effect, and equi- 


librium was quickly reestablished at differential salinity ~ — 0.01 ppm. 
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Fig. 13. Event #10 (Sep 1969): (a) dillerential MLD, (b) diflerential temperature, 


(c) differential salinity, and (d) wind speed 
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11. Event #11 (Oct 1969) 

In this period the MLD of case 2 was becoming deep again (h ~ 60m) and 
continued to deepen in this early winter regime (Fig. 7 on page 23). Increasing wind 
speeds and near Zero net heat values accounted for this effect. 

This return to the winter regime bidded a return to the differentia! MLD and 
differential temperature profiles seen in the first three events. Initial cooling of the 
shallower MLD of case 1 generated a negative spike in the differential temperature 
which moved positive as the decreased entrainment rate of case 2 inhibited entrainment 
of the layer. Further mixing acted to equalize MLD’s and return the differential MLD 
and differential temperature toward zero. After the differential temperature flux moved 
to a positive value it remained at a larger level than observed in the first three events, 
as it did in the summer regimes. This indicated that the mixed laver was in transition 
from the summer to wi-ter regimes and MLD's were shallower than those of the first 
cases. 

There was an immediate decrease in salinity in the first case, as in the other ex- 
periments, as the precipitation diluted the layer and this decrease continued until ces- 
sation of the precipitation event. Mixing by entrainment decreased the differential 
salinitv to a value of -0.2ppm, after which time it rapidly approached that of case 2. 

12. Event £12 (Nov 1969) 

This event was very similar in forcing conditions and resulting output to that 
of Event #11. In this event the MLD in case 2 was deeper than that of Event #11 as 
early winter deepening continued, and mean wind speeds appeared to be stronger with 
larger peak wind events (Fig. 7 on page 23). This deeper MLD vielded a larger differ- 
ential MLD after precipitation onset, as well as a lesser positive peak of differential 


temperature. This also yielded a simular differential salinity perturbation. 


C. EXPERIMENT SUMMARY 

The results of these events indicate that the effects of discrete precipitation events 
varied greatly as a function of low frequency seasonal heat flux forcing and higher fre- 
quency wind forcing. These two forcing boundary conditions strongly influenced the 
mixed laver thermal characteristics, as well as the salinitv. 

The average deeper, colder mixed laver in winter was significantly shallowed bv a 
single precipitation event. This shallowing allowed for the increased loss of the thermal 
energy of the shallower mixed laver to the colder atmosphere, supporting a temperature 


decrease in the mixed layer. Wind events that followed precipitation deepened the mixed 
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laver, but it did not deepen as much, nor as quickly as in the case without precipitation. 
Follow-on wind events deepened the mixed layer until entrainment returned the MLD 
to that of the case witkout precipitation. Entrainment had a similar effect on the dif- 
ferential mixed layer temperature. 

The effects of the precipitation event in winter were modulated by the discrete 
deepening and shallowing events that occurred just prior to precipitation onset. A 
deeper MLD increased the magnitudes of the differential values that occurred as a result 
of precipitation shallowing. A shallower MLD operated in the opposite sense to de- 
crease these values. Deepening and shallowing events occurring concurrent with pre- 
cipitation had a similar effect, acting to enhance or decreasing differential values. 

In summer the average shallower, warmer mixed layer was less significantly shal- 
lowed by these precipitation events. Given that the mixed layer was already relatively 
shallow, the decrease in MLD was not as significant as in winter. Significant thermal 
stratification in summer yielded a large increase in mixed laver temperature with precip- 
itation induced shallowing, as the thermal energy input at the top of the nuxed laver was 
concentrated into a shallower mixed laver. This increased mixed laver temperature 
produced a more stable mixed laver as the temperature jump at the base of the mixed 
laver increases. The increase in temperature jump tended to further inhibit entrainment 
and slowed deepening of the mixed laver. Deepening slowed enough that the differential 
temperature actually increased with time. 

The discrete precipitation events could have had a long term effect on the mixed 
laver. This was especially true in summer when MLD’s were shallower. In cases where 
strong winds generated strong mixing events, the mixed laver continued to demonstrate 
increased buoyancy, as long as 55 davs after the precipitation event, as the differential 


MLD became negative and differential temperature vielded a positive fluctuation. 
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V. PROCEDURES, RESULTS, AND ANALYSIS 


A. PROCEDURES 

The 400-day mixed layer simulation, as described previously, was conducted to ex- 
amine the accumulative effect of a time series of rainfall events in explaining salinity and 
temperature changes with time and depth. Discrete daily precipitation data was not 
available for this study. Only monthly totals were available. As a result a method was 
devised to simulate daily precipitation values using the available oceanographic and 
meteorological observations. With the lack of surface pressure and upper atmospheric 
data, precipitation prediction became difficult, thus prediction of localized atmospheric 
convective activity was difficult. Utilizing wind and cloud information, however, an at- 
tempt has been made to predict precipitation associated with frontal activity. This 
technique was based on meteorological principles developed by Petersen (1956). 

Specifically, this procedure utilized wind speed and direction, and cloud coverage 
data (Fig. 14 on page 42). With cloud cover data available in octaves, precipitation was 
allowed for cases with cloud coverage equal to or greater than 7/8. Considering pre- 
frontal and frontal rain, precipitation was allowed to occur for cases with wind directions 
from 000° (northerly Winds) clockwise to 220° (southwesterly winds), and wind speeds 
greater than 10 kts. This accounted for the southerly to southwesterly winds that pro- 
ceed a front and the “confused” northerly to southeasterly winds that accompany frontal 
passage. For post-frontal rain, precipitation was allowed to occur for cases with wind 
direction 220° (southwesterly winds) clockwise to 000° (northerly winds), and for 
Stronger wind speeds, greater than 25 kts. This accounted for the strong southwesterly 
to northwesterly winds that accompany a post-frontal regime and ensured consideration 
of all possible wind directions. Additionally, the post-frontal case required that the sea 
surface temperature be greater than the air temperature, as the cold frontal air passes 
over warmer water. Precipitation amount was adjusted to vield the 400-dav precipitation 
total of the observed monthlv data. The resulting rain with time is contained in Fig. 15 
on page 43. 

This “complex” alogarithm for precipitation forcing scheme was applied to the 
OPBL model at a range of possible amounts that may represent normal interannual 
variability in the North Pacific. A “medium” precipitation level was chosen bv adjusting 
total 400-day output to equal that of the summation of the 13 month period from 
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15. Complex Rain Model Output and Winds for 400 day period: (a) rain 
model output for 400 days from 23 Nov 1968 adjusted to equal actual 


monthlv totals, and (b) wind for same time period. 
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December 1968 to January 1970 (Table 3 on page 45). The “low” level was chosen to 
equal two-third's of the total 400-day precipitation total of the medium level case. The 
“high” level case was equal to one and one-half times the medium case and the “very 
high” case was made equal to three times the medium case (Fig. 16 on page 46). The 
same wind stress and heat flux atmospheric forcing values and initial conditions de- 


scribed in chapter II were used for all four model runs as in chapter IV. 


Table 2. PRECIPITATION TOTALS FOR THE 400-DA Y STUDY PERIOD 


quantitative 400-day total 
i Ra 








These four complex cases Were compared with the effects of “simple” constant 
evaporation and precipitation... Here the four “simple” precipitation and evaporation 
rates were specified to give constant evaporation and precipitation with the same total 
precipitation and evaporation as those of the four complex cases. These simple precip- 
itation cases were applied to the model so that the effects of discrete precipitation events 
could be compared to those of constant E-P. The resulting MLD, temperature, and 
salinity time series were differenced for comparison. The output mixed laver values of 
MLD, mixed laver temperature, and mixed laver salinity from complex and simple forc- 
ing were compared to the CTD observations. 

Advection was maintained in the model at 0.0 ppm/day. Considering the placement 
of advection m the saini equation: 


A S(E SP ASW 
ع تم ات‎ (5.1) 





It should be noted that an increase in advection would provide the same tendency for 


mixed layer salinity as would decreasing precipitation. 
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Table 3. MONTHLY PRECIPITATION AT OCEAN STATION "P" FOR 400 
DAY PERIOD 


Monthlv Precipitation for 400-day Period 





B. RESULTS AND ANALYSIS 
l. Complex Precipitation and Evaporation Forcing 

The response of the OPBL model depended upon the net precipitation for the 
cases of complex precipitation. The 400-dav total precipitation was increased from low 
to very high as indicated in Table 2 on page 44. An increase in total precipitation caused 
a general decrease in MLD's, as a result of increased buoyancy and decreased 
entrainment (Fig. 17 on page 47). As the precipitation level for complex precipitation 
forcing was increased from low to very high, it was larger in magnitude in winter than 
in summer, though the relative increase in summer was proportional to that of winter 
(Fig. 16 on page 46). This was due to the small amount of precipitation in summer 
compared to winter. This increase in precipitation in winter with increasing total pre- 
cipitation increased buoyancy more in winter than in summer, yielding shallower MLD's 
for larger complex precipitation amounts. Additionally, the largest magnitude of differ- 
ential MLD's occurred in winter (Fig. 17 on page 47 and Fig. 18 on page 48). 

Associated with the differentia] MLD flux that accompanied increased complex 


precipitation levels was a change in mixed layer temperature and salinity. Differential 
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Fig. 17. Complex Precipitation Forcing Predicted NILD’s: (a) solid line represents 
very high precipitation forced case and (b) dotted line represents low pre- 


cipitation forced case. 
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. 18. Subtracted Values of Complex Precipitation Forcing Output: All values 


are the result of subtracting Very Thigh precipitation forcing output from 
Low precipitation forcing, (a) diflerential MLD, (b) differential temper- 


ature, and (c) differential salinity. 
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Complex Precipitation Forcing Predicted Temperatures: (a) dotted line 
was low precipitation forced, (b) dashed line was medium precipitation 
forced, (c) chained-dashed was high precipitation forced, and (d) solid line 


was very high precipitation forced. 
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mixed laver temperature decreased onlv slightly in winter when the MLD was deep. 
Negative differential temperatures were induced bv the winter heat loss of the shallower 
mixed laver. As the summer approached, and MLD shallowed due to increased Q, and 
decreased winds, the differential MLD, had a more significant effect on differential 
temperature. Thus, the decreased MLD in the shallowing season concentrated more of 
the incident solar energy in a shallower layer than it would have without increased pre- 
cipitation. This produced a large positive differential temperature. 

Examination of the differential salinity showed a continual decrease in salinity 
as very high complex precipitation forcing inputed more fresh water into the mixed laver. 
Deepening events which yielded differential MLD’s produced differing entrainment 
events. Deeper MLD’s, characteristic of lower complex precipitation forced cases, en- 
trained water of higher salinity than in the higher complex case. Larger precipitation 
events input freshwater into the surface layer inhibiting entrainment, vielding a shallower 
MLD. This reduction in MLD further concentrated the fresh water, reducing surface 
S (Fig. 20 on page 51). 

From CTD observations, actual values of MLD, temperature, and salinity were 
compared to those values predicted by the OPBL model with complex precipitation 
forcing. Verification was hindered by data voids, especially during the first 95 days of 
the study. Nevertheless, it may be concluded that the MLD’s of the CTD data were best 
approximated by the model using very high complex precipitation. The envelopes of the 
fluctuaung MLD of all four complex precipitation forcing cases appeared to be deeper 
than those of the CTD casts. With the small number of CTD’s over any synoptic pe- 
riod, an hour-by-hour comparison was not possible, although MLD was predicted to 
fluctuate diurnally by as much as 50 m’day. As the level of precipitation was increased, 
the resulting MLD became shallower. Thus, the case with the largest precipitation 
amount appeared to most closely approximate the CTD envelope for MLD (Fig. 21 on 
page 52). 

The mixed laver temperatures compared favorably with observations in all four 
complex cases. As the study period progressed, the temperatures of the actual data in- 
creased more than those from complex precipitation forcing commensurate with the ac- 
tual MLD being shallower. Thus, the case with the highest complex precipitation 
produced the greatest shallowing, vielding the highest summer temperature. The peak 
temperature from the very high complex precipitation case, however, only slightly ex- 


ceeded that of the low precipitation case (Fig. 22 on page 53). 
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Fig. 21. Complex Precipitation Forcing Predicted MLD and CTD MLD: (a) the 
solid hne was very high precipitation forced and (b) the dotted line was 


low precipitation forced. 
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dotted line was low precipitation forced, (b) dashed line was mediuin pre- 
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Fig. 23. Complex Precipitation Forcing Predicted and CTD Salinity: (a) dotted 
line was low precipitation forced, (b) dashed line was medium precipi- 
tation forced, (c) chained-dashed was high precipitation forced, and (d) 


solid linc was very high precipitation forced. 
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The model-predicted salinity values vielded less certain results. Actual CTD 
salinity data behaved in a cyclic manner during the study period. Here the initial salinity 
of 32.5 ppm immediately increased with increasing MLD. As the observed MLD shal- 
lowed, the salinity decreased, until the minimum was achieved in summer. This was 
possibly associated with shallowing due to decreased wind and precipitation. As MLD 
deepened again around day 330, salinity increased (Fig. 23 on page 54). 

Overall, this observed cyclic behavior was best approximated by modeling with 
high complex precipitation forcing. The two cases with precipitation less than the high 
case yielded salinity values that consistently exceeded those of the actual data, and were 
non-cyclic. During winter deepening for these cases, the salinity increased dramatically 
from day 0 to day 100. During spring shallowing and into summer, the salinity contin- 
ued to increase gradually, and did not exhibited a cyclical decrease. For the very high 
complex case, the salinity exhibits a decreasing trend almost immediately. Here the 
buoyant damping of the turbulence by the strong precipitation generally countered the 
shear production of turbulence attributable to the wind stress. 

Though the high complex precipitation case did agree with the cyclic behavior 
of the CTD salinity data, there were salinity observations of significant disagreement 
with the model prediction. As stated previously, peak salinity in winter was less than the 
actual value. In summer the salinity minimum was larger than the observed minimum 
by ~ 0.7 ppm. 

2. Constant Precipitation and Evaporation Forcing 

The results obtained during this study were similar to those for the complex 
precipitation forcing study. Upon comparing the various magnitudes of simple precipi- 
tation forcing, it was expected that increasing levels of constant E-P would have a 
magnifying effect on the mixed laver parameters. For simplicity, high and low precipi- 
tation forcing cases were compared here. The increased level of precipitation in the high 
case yielded a shallower overall MLD than the low case (Fig. 24 on page 56). This 
shallowing effect was stronger in winter than in summer, as MLD’s tended to be deeper 
in winter. The peak predicted in the shallower period in the low case (Fig. 24 on page 
56 point a) did not appear in the high case. 

It appeared that winter temperatures were well verified for all simple and com- 
plex forced cases by the actual mixed laver temperatures, as the increased precipitation 


yielded only slightly more shallowing. In summer, however, differential temperatures 
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Simple Precipitation Forcing Predicted Temperature: (a) dotted line was 
low precipitation forced, (b) dashed line was medium precipitation forced, 
(c) chained-dashed was high precipitation forced, and (d) solid line was 


very high precipitation forced. 
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increased as the high case absorbed more solar heat in a shallower layer than the low 
case. 

As with the complex precipitation cases, the constant low and medium precipi- 
tation predictions vielded a proportional decrease in salinitv. The high precipitation case 
proved to be the most cyclic in time over the vear. The very high case yielded an almost 
monotonic decrease in salinity as a result of the large constant surface fresh water influx 
(Fig. 26 on page 59). 

From comparing simple to complex precipitation forcing, it was noted that the 
complex case produced relativelv more precipitation in winter than in summer, and there 
was greater evaporation in summer than in winter. The constant E-P forcing effect of 
the simple case thus had a larger precipitation rate during periods of no rain for the 
complex case (in winter when low winds occur and during the summer). In winter dif- 
ferential MLD fluctuated both positively and negatively about zero. It appeared to 
fluctuate negatively during a strong complex precipitation event yielding shallowing. It 
fluctuated positively during non-rain periods for the complex case, producing relative 
shallowing in the simple case. In winter these differential MLD values exceeded those 
of summer. 

Differential temperature in winter showed approximately equal temperatures for 
the four rain levels in both the simple and complex cases. Additionally, anv fluctuation 
in MLD had a relatively small effect on the differential temperature in the winter season. 
In summer a negative differential temperature was induced by precipitation being greater 
in the simple case than in the complex case. 

With respect to salinity, differential salinity tended to be positive in. winter. 
Here, the salinity of the complex case was usually less than that of the simple case due 
to the greater amount of precipitation in the complex case. It did fluctuate during this 
period, as periods of low winds vielded a decreased complex precipitation allowing sim- 
ple precipitation to exceed it and produce a negative fluctuation in differential salinity. 
In summer a positive d‘fferential salinity pre-dominated. Here, the lower precipitation 
for the complex case allowed the salinity of this case to increase. 

The very high simple precipitation case provided the best approximation to ac- 
tual MLD of all studv cases, both simple and complex. The large amount of precipi- 
tation produced consistently shallower MLD’s than the other simple cases. Additionally, 


it provided shallower MLD’s than occurred during the complex case when forced by low 
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Cases: All values are the result of subtracting complex high precipitation 
forcing output by Low precipitation forcing, (a) differential MLD, (b) 
difTerential temperature, (c) differential salinity, and (d) complex precipi- 


tation forcing. 
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winter winds and in all summer cases when low winds occur and MLD's tend to be 
shallowest (Fig. 28 on page 61). 

In winter, there was close agreement between all complex and simple forced 
cases and actual temperature. In summer, the shallowest MLD’s of the simple forced 
case produced the highest temperature, which though best approximated the actual 
temperature in summer, were still less than that of the actual case by 2.4°C. 

As with complex forcing, the high simple case provided the best approximation 
to the actual salinity values. The simple high rain case yielded a larger negative differ- 
ential salinitv in winter and a lower salinity in summer than the high complex case, while 
still retaining a cyclic behavior with a one-vear period. Its cvclic response was not as 


clear as in the complex case, however. 


VI. SUMMARY 


This study was conducted to examine the discrete precipitation event effects on the 
short-term and seasonal evolution of ocean mixed laver temperature and salinity struc- 
ture. Data for this study was for the years 1968 and 1969 at Ocean Station "P" located 


in the Northeast Pacific Ocean. Two kinds of experiments were performed: 


1. the first was to simulate the response of the mixed layer to discrete (isolated) pre- 
cipitation events, and 


2. the second was a quasi-realistic “complex” precipitation forcing experiment, re- 
presenting a realistic distribution in time of synoptic events over a 13-month pe- 
riod. 

The hypothetical experiment was conducted using twelve separate single discrete 
precipitation events over a 400-dav period. It examined the specific effects of a discrete 
precipitation event as a function of season without contamination bv other complex 
rainfall history. This experiment showed that the effect of a single discrete rain event 
can vary with season and the MLD (at onset of the rain event), as well as with the in- 
tensity of the wind stress. Additionally, it showed that a single rain event can have a 
lasting effect on mixed layer depth and temperature for up to 55 davs after the event 
Stops. 

For the main experiment a “complex” precipitation modeling scheme, producing 
discrete dailv precipitation events from surface meteorological observations was devised 
and tuned to vield realistic precipitation input for the mixed layer model. The precipi- 
tation model was adjusted to give differing amounts of total precipitation. Four different 
precipitation "intensities" were tested: 32 cm/400-day, 48 cm,400-day, 96 cm, 400-day, 
and 160 cin/400-day. The 50 cm’400-day case was chosen to yield the same net precip- 
itation as the observed precipitation data. These quasi-realistic precipitation forcings 
were applied to the mixed layer model. The resulting values of MLD, temperature, and 
salinity for this 400-day period were examined by comparing them with output values 
when the model was forced by “simple” constant net precipitation with the same total 
precipitation as in the complex cases. Finally, the values predicted by the model under 
Various precipitation forcing conditions were compared to.actual CTD values. 

In all study cases the MLD's approximated the MLD cycle of the actual data. The 


observed MLD’s were shallower than predicted by the model by 10m or more. Thus, the 
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model runs with the most precipitation best approximated the observed MLD's. It was 
noted that the greatest discrepancy occurred in the first winter season of the 400-day 
study period. However, a large data void existed during this time period. 

As with MLD, the temperatures of these studies generally approximated the tem- 
perature cycle and form of the actual data. The temperature in summer provided the 
greatest discrepancy in temperature magnitude with actual temperatures (exceeding 
those of the output data by at least ~ 2.5° C in the complex case and ~ 2.5°C in the 
simple case). 

Predicted salinity provided the greatest discrepancy between the modeled and the 
observed cycle. Rain rates of 32 cm/yr and 48 cm/vr (that equal to the observed total) 
vielded salinity values which were non-cvclic and, overall, higher in value than the ob- 
servations. The 96 cm'400-day case (both complex and simple) best approximated the 
salinity observations. Predicted salinity was less than observed in winter (by ~ 0.]0ppm 
in the complex case and ~ 0.lI ppm in the simple case) and greater than observed in 
summer (by ~ 0.07ppm in the complex case and ~ 0.07ppm in the simple case). The 160 
cm 400-dav precipitation cases produced a constant freshening of the mixed laver in 
both winter and summer. 

In the final analysis, while the discrete precipitation events had a different effect on 
the mixed laver than the simple cases, the difference was slight. Increased levels of pre- 
cipitation did enhance the difference between constant rain and the quasi-realistic rain 
cases, as seasonal differences in precipitation were amplified. The 96 cm/400-day 
quasi-realistic precipitation forcing case was selected as the most likely precipitation 
forcing for Ocean Station “P” during this time period. Despite magnitude discrepancies 
with observed rain at Station "P," the form and cycle of the MLD's and temperatures 
with respect to time were well approximated by this level of precipitation. The most 
important point for selection of this precipitation case was, however, the agreement with 
the actual salinity cycle. Here, the simple high precipitation case also approximated the 
shape of the salinity curve, but was not as cyclical, especially during the winter in 1969. 
This was probably because the complex precipitation case produced some periods with- 
out rain allowing for shallowing, and, therefore, decreased salinity. Even in winter, while 
the simple cases produced continuous rain. One such realistically modeled period acted 
to enhance cyclical behavior of salinity in the early winter of 1969, adding to MLD 


deepening and adding less fresh water during this period than that of the simple case. 


There are several reasons why observed mixed layer temperature and salinitv values 
may not agree with model predicted values: 
J. the inaccuracies and voids in CTD data, 


2. the inaccuracies in model forcing, especially in producing realistic complex precip- 
itation forcing, 


3. the inaccuracies in model initial conditions for T(z) and S(z), 
4. the inaccuracies in representing advection, and 


5. the mixed laver theory inaccuracies. 


Data voids do exist in the available CTD data at Ocean Station “P.” The largest 
data void occurred during the first 95 days of the study. At this time, there was an ob- 
servation at day 0 (Nov 23, 1968) and then no observations until the 95th day (Feb 26, 
1968). Additionally, data was sparse from days 95 to 180 and days 300 to 400. This 
would not be significant for data with small variance, as with the sea surface temper- 
ature. However, for data with strong variation such as the MLD and salinity data, the 
resulting analysis may be biased by lack of information at times of critical changes in 
mixed laver evolution. 

Inaccuracies In prec pitation forcing for the model mav exist. Four net precipitation 
rates, each multiples of the first, were chosen to bracket the most realistic precipitation 
amounts for Station "P.^ The fact that the 48 cm/400-dav case did not approximate the 
salinity cycle as well as the 96 cm/400-dav case, suggests that the observed precipitation 
rate may be lower than the actual amount. Allen (1963) outlined the significant number 
of potential rain error measurement sources onboard ship at Ocean Station “P.” These 
ere: 


]. the wind effects, 


t2 


the shadowing bv nearbv objects, 
the evaporation in the gauge, 


the sea spray entering the gauge, and 


^om سب‎ 


the tilting of the gauge orifice due to roll and pitch of the ship. 


Additionally, the rain model scheme used in this study was very crude. It may have 
inaccurately approximated the strength and duration of storm events, and would have 
been completely insensitive to local convection events (rain without wind). Such local 
convective events may account for a lower observed salinity, shallower MLD’s, and 


higher temperatures in summer. 
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Inaccuracies in meteorological data used in this study to compute heat and mo- 
mentum fluxes mav also exist. This data was collected in 1968 and 1969, and collection 
techniques have improved since. In addition to the wind stress and heat flux errors in- 
troduced directly into the model bv this data, these errors could produce unrealistic 
precipitation and evaporation events as well. 

Predicting horizontal advection is difficult, and given only surface meteorological 
data and scattered CTD data it 1s not possible with anv accuracy. Advective events have 
been known to occur at this station, as discussed in Chapter I. Such an event could 
greatly modify mixed laver structure. Advection was neglected for these studies, recog- 
nizing that an increase in advection of water of higher salinity has the same effect as 
decreasing rain and 1s necessarv to balance the salinity budget in the longer interannual 
term. This was another reason for examining various rain amounts. 

It is recommended that this experiment be repeated using observed discrete precipi- 
tation if possible. Alternatively, Tucker (1961) suggested a rainfall prediction scheme 
utilizing weather codes, using the standard international two-digit svstem. These 
follow-on studies could help answer the question as to the discrepancy in actual versus 
measured rainfall. Finallv, precipitation from such sources as satellite imaging could be 


interpreted as with Hakkarien and Adler (1987? and applied to the mixed laver model. 
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